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ABSTRACT 

The 3'-to-5' exoribonuclease in coronavirus (CoV) nonstructural protein 14 (nsp14-ExoN) 
mediates RNA proofreading during genome replication. ExoN catalytic residues are arranged in 
three motifs: I (DE), II (E), III (D). Alanine substitution of the motif I residues (AA-E-D, four 
nucleotide substitutions) in murine hepatitis virus (MHV) and SARS-CoV yields viable mutants 
with impaired replication and fitness, increased mutation rates, and attenuated virulence in vivo. 
Despite these impairments, MHV- and SARS-CoV ExoN motif I AA mutants (ExoN-AA) have 
not reverted at motif I in diverse in vitro and in vivo environments, suggesting that profound 
fitness barriers prevent motif I reversion. To test this hypothesis, we engineered MHV-ExoN-AA 
with 1, 2 or 3 nucleotide mutations along genetic pathways to AA-to-DE reversion. We show 
that engineered intermediate revertants were viable but had no increased replication or 
competitive fitness compared to MHV-ExoN-AA. In contrast, a low passage (P10) MHV-ExoN- 
AA showed increased replication and competitive fitness without reversion of ExoN-AA. 
Finally, engineered reversion of ExoN-AA to ExoN-DE in the presence of ExoN-AA passage- 
adaptive mutations resulted in significant fitness loss. These results demonstrate that while 
reversion is possible, at least one alternative adaptive pathway is more rapidly advantageous than 
intermediate revertants and may alter the genetic background to render reversion detrimental to 
fitness. Our results provide an evolutionary rationale for lack of ExoN-AA reversion, illuminate 
potential multi-protein replicase interactions and coevolution, and support future studies aimed at 


stabilizing attenuated CoV ExoN-AA mutants. 
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IMPORTANCE 

Coronaviruses encode an exoribonuclease (ExoN) that is important for viral replication, fitness, 
and virulence, yet coronaviruses with a defective ExoN (ExoN-AA) have not reverted under 
diverse experimental conditions. In this study, we identify multiple impediments to MHV-ExoN- 
AA reversion. We show that ExoN-AA reversion is possible but evolutionarily unfavorable. 
Instead, compensatory mutations outside of ExoN-AA motif I are more accessible and beneficial 
than partial reversion. We also show that coevolution between replicase proteins over long-term 
passage partially compensates for ExoN-AA motif I but renders the virus inhospitable to a 
reverted ExoN. Our results reveal the evolutionary basis for the genetic stability of ExoN- 
inactivating mutations, illuminate complex functional and evolutionary relationships between 
coronavirus replicase proteins, and identify potential mechanisms for stabilization of ExoN-AA 


coronavirus mutants. 
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INTRODUCTION 

The rapid evolution of RNA viruses represents a significant challenge for preventing, treating, 
and eradicating RNA viral diseases. High mutation rates in RNA viruses generate extensive 
opportunities to overcome evolutionary hurdles, such as antiviral drugs, host immunity, or 
engineered attenuating changes (1). The evolutionary pathways traversed by RNA viruses are 
shaped by natural selection, which will favor some evolutionary trajectories more than others 
based on whether mutations are beneficial, deleterious, or neutral (2). Predicting the likely results 
of RNA virus evolution is an important step for anticipating viral emergence and for developing 


escape-resistant antiviral drugs and vaccines (3, 4). 


Coronaviruses (CoVs) are a family of positive-sense RNA viruses that cause human illnesses 
ranging from the common cold to severe and lethal respiratory disease (5). All CoVs encode a 
proofreading exoribonuclease within nonstructural protein 14 (nsp14-ExoN) that is critical for 
replication, fidelity, fitness, and virulence, and ExoN-inactivation has been proposed as a 
strategy for live-attenuated vaccine development (6-15). As members of the DEDDh superfamily 
of exonucleases, CoV ExoNs hydrolyze nucleotides using four metal-coordinating amino acids 
arranged in 3 motifs: I (DE), II (E), III (D) (16, 17). Alanine substitution of ExoN motif I (DE- 
to-A A) disrupts ExoN biochemical activity in both SARS-CoV and human CoV 229E (hCoV- 
229E) (16, 18, 19). The betacoronaviruses murine hepatitis virus (MHV) and SARS-CoV 
tolerate disruption of ExoN activity [ExoN(-)] but display mutator phenotypes accompanied by 


defects in replication, competitive fitness, and evasion of innate immune responses (10, 13, 14). 
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EXON active site mutants in alphacoronaviruses, including transmissible gastroenteritis virus and 


hCoV-229E, have yet to be recovered and are proposed to be lethal for replication (19, 20). 


Given the critical role of ExoN in CoV biology and the elevated mutation rate, we expected that 
natural selection would repeatedly drive reversion of the ExoN-inactivating substitutions. In line 
with this expectation, ExoN motif III mutants of SARS-CoV and MHV rapidly and repeatedly 
revert ((14) and unpublished observations). In contrast, we have never detected partial or 
complete reversion of ExoN motif I mutants (ExoN-AA) in SARS-CoV or MHV during 10 years 
of study and hundreds of experiments. More specifically, we have not detected consensus or 
minority variants of any kind at the motif I AA codons in either virus strain during acute 
infections and prolonged passage in tissue culture and following treatment with multiple 
nucleoside analogues (6-11, 13, 14). SARS-CoV-ExoN-AA also is stable during acute and 
persistent animal infections in immunocompetent and immune-compromised mice (12). The lack 
of primary reversion is not due simply to reduced adaptive capacity, as both SARS-CoV- and 
MHV-ExoN-AA can adapt for increased replication (7, 14). Most strikingly, long-term passage 
of MHV-ExoN-AA (250 passages, P250) yielded a highly fit population that had directly 
compensated for defective proofreading through evolution of a likely high-fidelity RdRp (7). 
Yet, where primary reversion would have required just four total consensus mutations, MHV- 


ExoN-AA-P250 contained more than 170. 


In this study, we sought to determine whether specific genetic or fitness barriers prevent primary 


reversion of ExoN motif I AA. To this end, we identified and engineered viable genetic 
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pathways towards ExoN-AA motif I reversion in MHV (hereafter, ExoN-AA). Our results show 
that partial reversion did not confer a selective advantage compared to ExoN-AA. Further, 
ExoN-AA adapted within 10 passages for greater fitness than any of the intermediate revertants. 
Finally, restoration of WT-ExoN-DE in the setting of passage-selected mutations in the nsp12 
RNA-dependent RNA polymerase (RdRp) and nsp14-ExoN exacted profound fitness costs. 
Together, these data are the first observation of an ExoN(-) CoV genotype-fitness landscape and 
identify multiple genetic barriers underlying ExoN(-) motif I stability in MHV. Further, they 
suggest extensive coevolution between MHV replicase proteins during adaptation and reveal 


potential strategies for stabilizing ExoN mutant CoVs. 
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RESULTS 

Primary reversion of ExoN(-) motif I. MHV-ExoN(-), hereafter ExoN-AA, contains two 
engineered substitutions in each codon of motif I, such that complete reversion to WT-ExoN-DE 
requires mutations to all four sites (Figure 1A). Viral mutation rates in the absence of 
proofreading range from 10* to 10° mutations per nucleotide per round of replication (u) (1). 
Assuming an ExoN-AA mutation rate of 10 u and accounting for codon degeneracy, the 
probability of restoring the native amino acid sequence in a single round of replication is 108. 
Only rarely do ExoN-AA titers exceed 10° PFU/mL, so it is unlikely that ExoN-AA could 
navigate this genetic barrier in a single infectious cycle. Thus, we hypothesized that ExoN-AA 
reversion, if possible, would proceed incrementally. To identify potential pathways towards 
ExoN-AA reversion, we examined the possible single-nucleotide substitutions surrounding A89 
and A91 (Figure 1B). Three mutations are synonymous, and five mutations yield amino acids 
unlikely to coordinate with the positively-charged metals required for ExoN catalysis (glycine, 
valine, proline, threonine, and serine) (16, 19, 21, 22). One mutation per site can restore the 
acidic charge (i.e. AA-to-ED) but not the native amino acid. These variants have not been tested 
in a CoV EXON, but biochemical studies of E. coli DNA polymerase I ExoN mutants suggest that 
these conservative substitutions would not restore WT-like ExoN activity (23). We predicted 
stepwise pathways to ExoN-AA>DE reversion based on restoration of acidic charge followed 
by reversion to native amino acids (Figure 1C). We engineered and recovered variants in ExoN- 
AA requiring three mutations (3nt; ExoN-AD, ExoN-EA), two mutations (2nt; ExoN-DA, ExoN- 
ED, ExoN-AB), or one mutation (1nt; ExoN-DD, ExoN-EE) for reversion to WT-ExoN-DE 


(Table 1). We will hereafter refer to these mutants as intermediate revertants. All intermediate 
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revertants generated viable progeny during recovery, demonstrating that reversion to WT-ExoN- 
DE along these pathways is theoretically possible. The 3nt and 2nt mutants were genetically 
stable during recovery, as confirmed by dideoxy sequencing. However, both 1nt mutants (ExoN- 
DD and ExoN-EE) reverted to WT-ExoN-DE during three independent recovery attempts, 
suggesting that these two variants are less fit than WT-ExoN-DE and demonstrating that 
reversion by 1nt mutation is readily accessible. To test whether the 3nt or 2nt mutants would 
revert more rapidly than ExoN-AA (4nt), we passaged three lineages of each mutant 10 times at 
multiplicities of infection (MOI) of 0.5 and 0.01 PFU/cell. We harvested supernatants and 
screened for reversion by visual inspection of plaque phenotypes at each passage. WT-ExoN-DE 
and WT-like viruses produce uniform, large plaques, while ExoN-AA-like viruses yield small, 
variably-sized plaques (13). When we observed mixed plaque phenotypes, we sequenced three 
large plaques from each lineage to confirm reversion. The 3nt (ExoN-AD and ExoN-EA) and 2nt 
(ExoN-DA and ExoN-ED) intermediate revertants showed no evidence of reversion over 10 
passages at either MOI (Table 1). In contrast, the 2nt ExoN-AE contained WT-revertants by P2 
in all lineages at MOI = 0.5 PFU/cell and by P8 in one lineage at MOI = 0.01 PFU/cell. Once 
observed, WT-revertants dominated the ExoN-AE population for the remaining passages. These 
data indicate that at least one 2nt mutation pathway can lead to full reversion in tissue culture. 
The probability of ExoN-AE arising during a single infectious cycle of ExoN-AA is low but 
theoretically achievable (-10?), so ExoN-AA could conceivably revert within just two infectious 
cycles. However, complete reversion has never been observed even during prolonged passage or 
persistent infections, suggesting that additional barriers to the replication, fitness, or maintenance 


of intermediate revertants exist. 
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Partial reversion of MHV-ExoN(-) motif I does not confer a selective advantage. Because 
the intermediate revertants are viable as recombinants but are not found in ExoN-AA 
populations, we hypothesized that they confer no selective advantage over ExoN-AA (8, 9, 13). 
To test this hypothesis, we first analyzed replication of the 3nt and 2nt intermediate revertants 
(Figure 2A). All variants achieved similar peak titers to ExoN-AA, but detailed examination of 
their kinetics suggested a potential delay of up to 1.5 hours for all intermediate revertants 
compared to ExoN-AA. Of note, ExoN-AE was the most delayed, and we detected WT-ExoN- 
DE revertants in two of three replicates, suggesting increased selective pressure against this 
variant. We next measured the competitive fitness of each intermediate revertant relative to a 
recombinant ExoN-AA containing seven silent mutations in the nsp2 coding region (ExoN-AA- 
reference). Intermediate revertants were mixed with an equal titer of ExoN-AA-reference at a 
combined MOI = 0.05 PFU/cell and passaged four times. The ratio of each intermediate revertant 
to ExoN-AA-reference was quantified at each passage by RT-qPCR, and the change in ratio over 
time was used to calculate their relative fitness. WT-ExoN-DE was significantly more fit than 
ExoN-AA, whereas the intermediate revertants (ExoN-AD, -EA, -DA, and -ED) had no 
increased fitness relative to ExoN-AA (Figure 2B). The apparent increased fitness of ExoN-AE 
resulted from all lineages reverting to WT-ExoN-DE during the experiment. Finally, our 
previous studies have shown that adaptation of ExoN-AA includes partial compensation for the 
replication fidelity defect, as measured by reduced susceptibility to the mutagen 5-fluorouracil 
(5-FU) (7-11, 24). None of the intermediate variants demonstrated statistically significant 
differences in 5-FU sensitivity as compared to ExoN-AA (Figure 2C). Thus, with the exception 


of the ExoN-AE>DE revertants, no 3nt and 2nt intermediate genotypes along our predicted 
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pathway demonstrated an advantage in replication, fitness, or fidelity that would favor their 
maintenance or expansion in the viral population. Thus, natural selection is unlikely to drive 


ExoN-AA down these pathways towards reversion. 


Secondary adaptations outside of ExoN-AA motif I increase fitness along alternative 
pathways. Although we did not find fitness advantages to intermediate revertants, we also did 
not identify profound fitness costs that would drive their immediate loss from populations. We 
have previously demonstrated that during 250 passages (P250), ExoN-AA can adapt for 
increased replication, fitness, and fidelity via secondary mutations outside of motif I (7). We 
tested whether secondary adaptive mutations could exceed the fitness of ExoN-AA and its 
intermediate revertants. To examine the early adaptation of ExoN-AA, we studied passage 10 
from the P250 passage series (Figure 3). ExoN-AA-P10 retains the ExoN-AA motif I genotype 
but has increased replication and reduced susceptibility to 5-FU, altogether manifesting in 
greater relative fitness (Figure 3) (7). We identified only six total mutations within ExoN-AA- 
P10 by dideoxy sequencing (Table 2), indicating that rapid adaptation of and compensation for 
ExoN-AA requires relatively few genetic changes at the consensus level. To test whether 
interactions between multiple mutations or population level effects contribute to ExoN-AA-P10 
fitness, we isolated a plaque-purified clone of ExoN-AA-P10. The clone replicated to higher 
titers than the ExoN-AA-P10 population but had identical 5-FU sensitivity and relative fitness 


(Figure 3), indicating that genomes derived from a single virus plaque encode the adaptive 


changes required by the total population. Together, these data demonstrate that mutations outside 


of ExoN(-) motif I can confer greater fitness advantages than intermediate revertants even at 


10 
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early passages. These early adaptive mutations likely reduce the selective pressure for motif I 


reversion and place the intermediate revertants at a selective disadvantage. 


Adaptive mutations in nsp12 and nsp14 that increase ExoN-AA fitness confer significant 
fitness costs to WT-ExoN-DE. Mutational fitness effects are highly dependent upon the genetic 
background (25-27). In addition to reducing selective pressure for reversion, mutations 
conferring increased fitness to ExoN-AA might also reduce the benefits of motif I reversion. We 
previously reported that long-term passage of ExoN-AA selects for secondary adaptive 
mutations in the nsp12 RdRp and nsp14 (nsp12-P250 and nsp14-P250) (7). Nsp12-P250 contains 
7 nonsynonymous mutations that partially compensate for defective proofreading and increase 
ExoN-AA fitness. Nsp14-P250 contains 6 nonsynonymous mutations, including a conservative 
D-to-E substitution in ExoN motif III, and increases ExoN-AA fitness without compensating for 
defective proofreading. To test whether the fitness effects of passage-associated mutations in 
nsp12-P250 and nsp14-P250 depend upon the ExoN-AA genotype, we engineered a WT motif I 
(ExoN-DE) into viruses containing nsp12-P250 and nsp14-P250, alone and together, and 
analyzed replication, 5-FU sensitivity, and competitive fitness. Compared to WT-ExoN-DE, both 
ExoN-DE-nsp12-P250 and ExoN-DE-nsp14-P250 displayed delayed and decreased replication 
(Figure 4A). In 5-FU sensitivity assays, ExoN-DE-nsp14-P250 was indistinguishable from WT- 
ExoN-DE, while both variants containing nsp12-P250 (ExoN-DE-nsp12-P250 and ExoN-DE- 
nsp12/14-P250) were significantly more sensitive to 5-FU (Figure 4B). Finally, the nsp12-P250 
and nsp14-P250 mutations significantly decreased fitness relative to WT-ExoN-DE (Figure 4C). 


We detected no statistical differences between the specific infectivity of WT-ExoN-DE and any 
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of the nsp12-P250 and nsp14-P250 variants in isolated infections (Figure 4D). Thus, mutations 
in nsp12 and nsp14 that arose in the ExoN-AA background were detrimental to replication, 
mutagen sensitivity, and competitive fitness in the presence of a fully-reverted ExoN-DE. These 
results support the conclusion that the adaptive pathways available to ExoN-AA may stabilize 
the ExoN-AA genotype, reducing both the selective pressure for, and the potential benefits of, 


primary reversion. 
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DISCUSSION 

In this study, we demonstrate that the stability of the ExoN(-) motif I genotype in MHV (ExoN- 
AA) is a consequence of the limitations and opportunities of the genetic landscape it explores 
during replication (Figure 5). Our results support a model in which the viable adaptive pathways 
leading to direct reversion of motif I from AA-to-DE are relatively flat on a fitness landscape, 
with intermediate revertants remaining phenotypically ExoN(-) and conferring no fitness 
advantage over ExoN-AA. In contrast, at least one alternative adaptive pathway is readily 
accessible and imparts immediate fitness gains over ExoN-AA. We propose that even minimal 
alternative pathway adaptive fitness gains reduce the likelihood and benefits of motif I reversion, 
until eventually the changing genetic background renders reversion detrimental. These data and 
this model suggest that selection during replication favors immediate, incremental fitness gains 
along the most accessible pathway rather than dramatic fitness increases across a larger genetic 
barrier. While this study focused on the issue of primary ExoN-AA reversion, it raises intriguing 
questions about the remaining topography of the ExoN-AA fitness landscape. Ongoing studies in 
our laboratory will determine whether ExoN-AA-P10 (and P250) represent a conserved pathway 
to compensation for defective proofreading or if CoV genomes accommodate multiple solutions 
to ExoN(-)-associated debilitations. Supporting the possibility of limited potential pathways, a 
recent study of vaccine-derived polioviruses determined that evolution for virulence in parallel 
epidemics proceeded through a limited set of genetic changes rather than a plethora of distinct 


trajectories (4). 
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A key finding of this study is that mutations compensating for ExoN-AA are detrimental to WT- 
ExoN-DE, yet the mechanisms underlying this disparity are not clear. We expected that 
compensatory mutations in nsp14-250 would hobble WT-ExoN-DE by inducing structural 
changes in the active site. However, the replication and fitness defects observed in ExoN-DE 
viruses containing nsp14-250 were not associated with altered proofreading, as measured by 5- 
FU sensitivity (Figure 4). In fact, the nsp14-250 mutations do not substantially affect 5-FU 
sensitivity in either the ExoN-AA or WT-ExoN-DE backgrounds (7), suggesting that these 
mutations enhance some other stage of the MHV replication cycle, perhaps through interactions 
with other replicase proteins or with viral RNA. If so, these results imply that the proteins of the 
MHYV replication holoenzyme must evolve in a highly cooperative manner to compensate for 
defective proofreading while also maintaining their intricate and intimate functional 
relationships. In line with this hypothesis, the likely high-fidelity nsp12-250 disrupted 
replication, 5-FU resistance, and fitness in the presence of an intact ExoN-DE catalytic motif. 
Given the inverse relationship between polymerase fidelity and replication speed (the slower the 
polymerase proceeds, the more time it has to discriminate incoming nucleotides) (28-31), nsp12- 
250 could fracture the delicate kinetic balance between polymerization and nsp14-ExoN- 
mediated excision during replication. Support for this hypothesis comes from a recent study 
establishing that the SARS-CoV-nsp12 replicates with lower fidelity than the RdRp of dengue 
virus and from the observation that the likely high-fidelity RdRp variant MHV-nsp12-V553I 
delays replication of proofreading-capable MHV (9, 32). As with nsp14-250, nsp12-250 
mutations could have additional effects on RdRp functions or protein-protein interactions that are 


incompatible with an active ExoN domain. Finally, it is possible that the adaptive mutations in 
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nsp12-250 and nsp14-250, including synonymous changes, could modify RNA sequence features 
or secondary structures in ways that interfere with WT-ExoN-DE. Defining the molecular 
interactions underlying the strain-specificity of ExoN-AA-compensatory mutations using 
recombinant viruses or reconstituted biochemical systems will provide a unique opportunity to 
explore the complex evolutionary and functional relationships that underpin the coronavirus 


replication machinery. 


Our results also extend existing studies of CoV ExoN motif I. Motif I AA>DE mutations in the 
SARS-CoV nsp14-ExoN dramatically reduce nuclease activity in biochemical assays, but no 
study has examined the contributions of each residue independently (16, 18). While we cannot 
exclude the possibility that intermediate revertants of ExoN-AA retain some level of nucleolytic 
activity, the lack of consistent or statistical differences in replication, 5-FU sensitivity, or 
competitive fitness relative to ExoN-AA supports previous studies demonstrating that motif I DE 
is essential for ExoN function (16, 18, 19). Given these results, we were surprised to observe 
repeated reversion of the ExoN-AE but not the other two 2nt variants, ExoN-DA and ExoN-ED. 
One potential explanation is that the specific mutational bias of ExoN-AE makes the revertant 
mutations more accessible than in ExoN-DA or ExoN-ED. Alternatively, if ExoN-AE has 
profound replication or fitness defects, selection could drive primary reversion more quickly 
away from this genotype. Consistent with this hypothesis, ExoN-AE reverted more quickly at a 
higher MOI, where natural selection acts more efficiently on a larger population size (Table 1) 


(33). Biochemical studies of the ExoN-AA intermediate revertants will be valuable to determine 
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whether nsp14-ExoN-AE differs from the other intermediate revertants and to define the 


structural and kinetic features of nsp14-ExoN catalysis. 


Finally, our studies suggest that compensatory mutations identified through long-term passage 
could be used to stabilize the ExoN-AA genotype. In particular, the high-fidelity nsp12-P250 
could reduce the probability of reversion by reducing mutational sampling within motif I (34), 
and both nsp12-P250 and nsp14-P250 render the MHV genome inhospitable to a WT-ExoN-DE. 
Together, these studies argue that experimental evolution can generate reagents to define critical 
interactions involved in CoV replication and can identify new strategies for stabilizing attenuated 


CoVs. 
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MATERIALS AND METHODS 

Cell culture. Delayed brain tumor (DBT-9) cells (35) and baby hamster kidney 21 cells 
expressing the MHV receptor (BHK-R) (36) were maintained at 37°C in Dulbecco’s Modified 
Eagle Medium (DMEM, Gibco) supplemented with 10% serum (HyClone FetalClone II, GE 
Healthcare or Fetal Bovine Serum, Invitrogen), 100 U/mL penicillin and streptomycin (Gibco), 
and 0.25 uM amphotericin B (Corning). BHK-R cells were further supplemented with 0.8 
mg/mL G418 selection antibiotic (Gibco). The infectious clone of the murine hepatitis virus 
strain A59 (MHV-A59; GenBank accession number AY910861) was used as the template for all 


recombinant viruses . 


Determination of viral titer by plaque assay. Virus samples were serially diluted and 

inoculated on subconfluent DBT-9 cell monolayers in either 6- or 12-well format. Cells were 
overlaid with 1% agar in DMEM and incubated overnight at 37?C. Plates were fixed with 4% 
formaldehyde and agar plugs were removed. The number of plaques per well was counted by 


hand and used to calculate titer (36). 


Plaque purification of viral populations. DBT cells were infected with serial dilutions of virus 
and overlaid with 1% agar in DMEM. Single plaques were isolated with glass Pasteur pipettes, 
resuspended in PBS containing calcium and magnesium, and inoculated onto fresh DBTs. This 


process was completed 3 times before generating experimental stocks. 
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Cloning and recovery of recombinant viruses. Site-directed mutagenesis in MHV genome 
fragments was performed using “round the horn” PCR (originally described in (37)). Briefly, 
adjacent primers containing the mutation of interest were 5'-phosphorylated using T4 
polynucleotide kinase (NEB, M0201S) using the buffer from the T4 DNA ligase, which contains 
ATP (M0202S). PCR was performed on a plasmid template using the Q5 High-fidelity 2x 
Master Mix (NEB, M0492L), with primers at final concentration of 500nM. The linear 
amplification product was purified using the Promega Wizard SV Gel and PCR Clean-up System 
(Promega Corporation, A9282), and 4 uL was ligated at 16°C overnight with the T4 DNA ligase 
(NEB M02025). After transformation into chemically-competent Top10 E. coli (lab-derived) and 
expansion in liquid culture, the MHV segment of each plasmid was sequenced. Viruses were 
constructed, rescued, and sequenced as described previously (7, 13, 36). Experimental stocks 
were generated by infecting a subconfluent 150 cm? flask of DBT-9 cells at MOI of 0.01 
PFU/cell. Flasks were frozen at -80°C when monolayers were fully involved, approximately 20- 
28 hours post-infection depending on the variant. After thawing, the supernatant was clarified by 
centrifugation at 4,000 x g (Sorvall RC 3B Plus; HA-6000A rotor) for 10 min at 4°C. For 
intermediate revertants, stocks were generated in serum-free DMEM and processed as above 
before being concentrated roughly 10-fold by centrifugation at 4,000 x g using Amicon Ultra-15 
Centrifugal Filter Units, 100kDa (EMD Millipore, UFC910008). The virus titer of each stock 


was determined by plaque assay using DBT-9 cells as described above. 


Passage of ExoN intermediate revertants. Intermediate revertants of ExoN-AA were passaged 


10 times on subconfluent DBT-9 cell monolayers in 24-well plates at an estimated MOI of either 
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0.01 or 0.5 PFU/cell. Supernatants were harvested at 24 and 20 hours post-infection for MOI = 
0.01 and 0.5 PFU/cell, respectively, and screened for WT reversion by plaque assay. At least 


three WT-like plaques were sequenced for each lineage to confirm motif I reversion. 


Replication kinetics. Viral replication kinetics in DBT-9 cells were determined at indicated 
MOIs as described previously (11). Replicates were synchronized by 30-minute incubation at 
4°C before transferring to the 37°C incubator. Supernatant (300 uL) was harvested at the 


indicated time points and titered by plaque assay. 


Determination of specific infectivity. Subconfluent monolayers of DBT-9 cells in 24-well 
plates were infected with the indicated virus at MOI = 0.05 PFU/cell, and supernatant was 
harvested at 16 hours post-infection. Genomic RNA in supernatant was quantified using one-step 
reverse transcription quantitative RT-PCR (RT-qPCR) on TRIzol-extracted RNA as described 
previously (9). Briefly, genomic RNA was detected with a 5’ 6-carboxyfluorescein (FAM) and 
3” black hole quencher 1 (BHQ-1) labeled probe targeting nsp2 (Biosearch Technologies, 
Petaluma, CA), and RNA copy number was calculated by reference to an RNA standard derived 
from the MHV A fragment. Samples were plated in technical duplicate to minimize well-to-well 
variation. Titers were determined by plaque assay in DBT-9 cells, and specific infectivity was 


calculated as PFU per supernatant genomic RNA copy. 


5-fluorouracil sensitivity assays. Stock solutions of 5-fluorouracil (Sigma F6627) were 


prepared in dimethyl sulfoxide (DMSO). Sensitivity assays were performed in 24-well plates at 
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MOI = 0.01 PFU/cell, as previously described (7). Cells were incubated with drug for 30 minutes 
prior to infection. Supernatants were harvested at 24 hours post-infection, and titers were 


determined by plaque assay. 


Competitive fitness assays. ExoN-AA-reference and WT-ExoN-DE-reference viruses were 
marked with 7 consecutive silent mutations within nsp2 (wild-type: 1301-TTCGTCC- 1307; 
reference: 1301-CAGCAGC-1307) by round the horn PCR, as described above. Competitions 
were performed in triplicate on DBT-9 cells in 12-well plates, plated at a density of 1 x 10° cells 
per well 24 hours prior to infection. Cells were infected at a total MOI of 0.1 PFU/cell (MOI = 
0.05 PFU/cell each for competitor and reference virus). Supernatants were harvested 15 and 16 
hours post-infection for experiments with ExoN-AA-reference and WT-ExoN-DE-reference, 
respectively, and passaged 4 times. Samples were titered between all passages to maintain total 
MOI of 0.1 PFU/cell. RNA was extracted from 70 uL of supernatant using QIAamp 96 virus 
OlAcube HT kit on the QIAcube HT System (Qiagen). Each RNA sample was analyzed by one- 
step RT-qPCR with two SYBR Green assays. Reference viruses were detected with forward 
primer SS-qPCR-Sil-F (5'-CTATGCTGTATACGGACAGCAGT-3"; 200nM final) and reverse 
primer SS-qPCR-R2 (5'-GGTGTCACCACAACAATCCAC-3', 200nM final). Competitors were 
detected with forward primer SS-qPCR-WT-F (5'-CTATGCT-GTATACGGATTCGTCC-3", 
450 nM final) and reverse primer SS-qPCR-R2 (5'-GGTGTCAC-CACAACAATCCAC-3', 450 
nM final). RNA samples were diluted 1:100 prior to RT-qPCR with Power SYBR Green RNA- 
to-Ct 1-step kit (Applied Biosystems) according to the manufacturer's protocol. Duplicate wells 


were averaged, and values were excluded from subsequent analysis if the duplicate wells differed 


20 


1san6 Aq eLOz “pz Amp uo /B10 use 1A(//: dy wo pepeojuwoq 


(0) 
E 
G 
O 
(05) 
iv 
U 
O 
jae 
D 
(o 
"te 
O 
N 
=) 
[en 
= 
me) 
E 
t. 
(15) 
O 
O 
< 


Journal of Virology 


Journal of Virology 


389 


390 


391 


392 


393 


394 


395 


396 


397 


398 


399 


400 


401 


by > 0.5 Ct. The relative abundance of competitor and reference were determined by subtracting 
Ct thresholds (ACteompetitor = Ctcompetitor ^ Ctreference) and converted to reflect the fold-change in 


940! competitor) "The log ¡yAratio was plotted against passage number, and the change 


ratio (Aratio = 
in log ¡pAratio (i.e. slope of linear regression) is the relative fitness. Note that regressions were fit 


only through P1-P4, as slight deviations in 1:1 ratio in the input (PO) can skew the slope. 


Statistical analysis. GraphPad Prism 6 (La Jolla, CA) was used to perform statistical tests. Only 
the comparisons shown [e.g. ns or asterisk(s)] within the figure or legend were performed. In 
many cases the data were normalized to untreated controls. This was performed using GraphPad 


Prism 6. The number of replicate samples is denoted within each figure legend. 
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FIGURE LEGENDS 


Figure 1. Sequence landscape around ExoN-AA motif I. (A) ExoN motif I nucleotide sequences. (B) 
Landscape of single-nucleotide substitutions within ExoN-AA motif I. (C) Predicted pathways to 
reversion of ExoN-AA. Variants marked with # reverted to WT during three independent recovery 


attempts. 


Figure 2. Intermediate revertants of ExoN-AA motif I do not have selective advantages. (A) 
Replication kinetics at MOI = 0.01 PFU/cell plotted as mean + SD of n = 3. (B) Competitive fitness of 
each variant relative to ExoN-AA. Viruses were competed with a tagged ExoN-AA-reference strain, and 
relative fitness was normalized to the mean of ExoN-AA. (C) 5-fluorouracil sensitivity at MOI = 0.01 
PFU/cell. Statistical significance of each variant relative to ExoN-AA was determined by one-way 
ANOVA with multiple comparisons (Panel D) two-way ANOVA with Dunnett’s multiple comparisons 
(panel C). ****p<0.0001; ns = not significant. Data in (B) and (C) represent mean + SD of n = 6. Boxed 
values have the same significance. "All lineages of ExoN-AE reverted to WT-ExoN-DE during the 


experiment. 


Figure 3. ExoN-AA adapts for increased fitness within 10 passages. (A) Replication kinetics of 
indicated viruses at MOI = 0.01 PFU/cell plotted as mean + SD of n = 3. (B) 5-fluorouracil sensitivity at 
MOI = 0.01 PFU/cell. (C) Competitive fitness of individual recombinants relative to ExoN-AA. Viruses 
were competed with a tagged ExoN-A A-reference strain, and relative fitness was normalized to the mean 
of ExoN-AA. Statistical significance of each virus relative to EXoN-AA was determined by two-way 


ANOVA with Dunnett's multiple comparisons (panel B) or by one-way ANOVA with multiple 
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comparisons (Panel C). ****p < 0.0001, ns = not significant. LOD = limit of detection. Data in (B) and 


(C) represent mean + SD of n = 6. Boxed values have the same significance. 


Figure 4. Mutations that increase ExoN-AA fitness are detrimental in the presence of WT-ExoN- 
DE. (A) Replication kinetics of indicated viruses at MOI = 0.01 PFU/cell plotted as mean + SD of n =3. 
(B) 5-fluorouracil sensitivity at MOI = 0.01 PFU/cell, mean + SD of n = 6. (C) Competitive fitness of 
individual recombinants relative to WT-ExoN-DE. Viruses were competed with a tagged WT-ExoN-DE 
reference strain, and relative fitness was normalized to the mean of WT-ExoN-DE, mean + SD of n = 6. 
(D) Specific infectivity (genomes per PFU) from isolated infections, mean + SD of n = 4.. Statistical 
significance of each virus relative to WT-ExoN-DE was determined with two-way ANOVA with 
Dunnett's multiple comparisons test (panel B) or by ordinary one-way ANOVA with Dunnett's multiple 


comparisons test (panels C and D). **p < 0.01; ***p < 0.001; ****p « 0.0001; ns = not significant. 


Figure 5. Model for the in vitro evolution of MHV-ExoN-AA. MHV-ExoN-AA (black dot) is a low- 
fitness variant. Reversion to WT-ExoN-DE would dramatically increase fitness but can only be achieved 
by traversing a flat landscape and climbing a steep fitness cliff (dotted white arrows). However, 
secondary mutations that incrementally increase fitness are more accessible (solid white arrow). 
Eventually, the genetic background changes enough that reversion becomes detrimental (dotted yellow 


line). 
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TABLE FOOTNOTES 
Table 1. Recovery and passage of intermediate revertants. 


n.d.: not done. 
“Bolded nucleotides must mutate to reach WT-ExoN-DE genotype. 
Recovered viruses were subjected to 10 passages at the indicated MOI. Samples were screened for wild- 


type revertants by plaque assay, and revertant lineages were sequence-confirmed. 


Table 2. Mutations in ExoN(-) P10. Data derived from dideoxy sequencing. 


“Mutation present at approximately 50% of population. 


ÞMHV HE is not transcribed in tissue culture. 


“Amino acid numbers designate positions within cleaved nsps, not the polyprotein. 
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Table 1. Recovery and passage of intermediate revertants. 

n.d.: not done. 

“Bolded nucleotides must mutate to reach WT-ExoN-DE genotype. 

Recovered viruses were subjected to 10 passages at the indicated MOI. Samples were screened for wild- 


type revertants by plaque assay, and revertant lineages were sequence-confirmed. 


# of reverted lineages by 
passage 10° 


Accepted Manuscript Posted Online 


# of mutations to Motif I —_ 

Virus WT-ExoN-DE sequence? MOI = 0.01 MOI = 0.5 
ExoN-AA 4 GCA...GCT 0/3 0/3 
ExoN-AD 3 GCA...GAT 0/3 0/3 
ExoN-EA 3 GAA...GCT 0/3 0/3 
ExoN-DA 2 GAT...GCT 0/3 0/3 

ExoN-AE 2 GCA...GAA 1/3 (by P8) 3/3 (by P2) 
9 ExoN-ED 2 GAA...GAT 0/3 0/3 
E ExoN-EE 1 GAA...GAA n.d. n.d. 
3 ExoN-DD 1 GAT...GAT n.d. n.d. 

= 

3 WT-ExoN-DE 0 GAT...GAA n.d. n.d. 
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Table 2. Mutations in ExoN(-) P10. Data derived from dideoxy sequencing. 
“Mutation present at approximately 50% of population. 
"MHV HE is not transcribed in tissue culture. 


“Amino acid numbers designate positions within cleaved nsps, not the polyprotein. 
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Mutation Nucleotide Amino acid 
number change Protein change” 
1 G2520A* nsp2 D524N 


A3080G*  nsp3 Silent 
T16017A  nspl2  M814K 
A17836G^ nspl3 1492M 
G22673A ^ HE” noncoding 
A29208CC M Silent 
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